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Numerical Solution of Supersonic Laminar Flow over 
an Inclined Body of Revolution 


C. M. Hung* 

NASA Ames Research Center, Moffett Field, Calif 


A mixed explicit-implicit scheme is used lo solve the time-dependent thin-layer approximation of the Navier- 
Stokes equations for a supersonic laminar flow over an inclined body of revolution. Test cases for Mach 2.8 flow 
over a cylinder with 15-deg flare angle at angles of attack of 0, 1, and 4 deg are calculated. Good agreement is 
obtained between the present computed results and experimental measurements of surface pressure. A pair of 
vortices on the leeward and a peak in the normal force distribution near the flared juncture are predicted; the 
role of circumferential communication is discussed. 


Nomenclature 


A 

= right-hand side of Eq. (5) 

Cfx 

= skin-friction coefficient along 
streamwise direction 


= specific heat at constant volume 

D 

= diameter of cylindrical portion 

E 

= total energy, -1-0.5 

FGM 

= transport inviscid flux vectors, Eq. 
(I) 


= transport viscous flux vector, Eq. (1) 

L 

= length of cylindrical portion 


= split finite difference operators in 


= hyperbolic and parabolic operators in 
Tf direction 

A/« 

= freestream Mach number 

Pr 

= Prandtl number, 0.72 

P 

= static pressure, (7- l)pe, 

P, 

= freestream static pressure 

Q 

= forcing terms in Eq. (1) 

QiiQi*Qv 

= split forcing terms in 

R 

= gas constant 

Re, 

= Reynolds number based on length L 

r 

= radial direction 

rb(x),r^{x),r^(x) 

= radius of body, fine mesh, and outer 
coarse mesh, respectively 

t 

= time 


= time increment of outer and inner 
meshes 

U 

= conservative variables (p,pw,pu, 
pw,pE) 

U' 

= variables in operator {u\v' , 


= freestream velocity 

u,v,w 

= velocity components in (x,r,<t>) 
coordinates 

u',v' 

= transformed velocity components, 
w ' = t; + viand v' = uianO, 

X 

= axial direction 
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= angle of attack 
= ratio of specific heat, 7 = 1 .4 
= minimum fine mesh spacing 
= flare angle 

=tan - * ) local body slope angle 

= second viscosity coefficient, X = 

-Vifi 

- molecular viscosity coefficient 
= viscosity parameter in Eq. (5) 

= transformed cylindrical coordinates 
= density 

= thin-layer approximated viscous 
terms 

= energy dissipation and diffusion 
terms 

Introduction 

S upersonic now over a body of revolution at angle of 
attack is of interest to both research and design engineers. 
A typical problem is to predict the nature of the three- 
dimensional flowfield around a cylinder with flared afterbody 
(Fig. I). This configuration is seen in missiles, rockets, and 
space launch vehicles, in which the flare provides a natural 
fairing for multistage spacecraft or a simple means for 
shifting the center of pressure rearward and thereby im- 
proving the margin of static stability. It can also be used to 
simulate the underexpansion plume of rockets. 

The pressure rises, due to the flare-shock, induce flow 
recirculation near the corner. At zero angle of attack, the flow 
is two-dimensional (axi symmetric) and the recirculation is 
closed as a toroid, provided the flow is steady. At angle of 
attack, the flare shock is stronger on the windward and 
weaker on the leeward. This variation of shock strength 
results in a pressure variation in the circumferential direction 
which induces a strong crossflow and makes the shock-wave 
and boundary-layer interaction highly three-dimensional. As 
the angle of attack increases, the crossflow first separates 
inside the recirculation region near the corner of the lee side. 
As the angle of attack is increased further, the crossflow also 
separates around the cylinder and the flare. Figure 2 shows a 
china clay laminar flow pattern and a sketched three- 
dimensional separation model experimentally observed by 
Ericsson et al.* and Fig. 3 shows a fluorescent oil 
visualization of a turbulent flow on a 23 deg cone frustum by 
Chyu and Coe.^ Both figures Indicate a circumferential 
communication and a pair or pairs of vortices shedding into 
the lee-side "‘wake. 

Because of the complexity, previous theoretical and ex- 
perimental studies considered the flare and Ice-sidc 
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separations separately. The flare recirculation is studied as 
flow over a compression corner and the lee-side separation 
is modeled as a flow over a cone.®*® To understand the 
detailed structure of the complete flowfleld, however, the 
problem should be considered as a whole, and a complete 
three-dimensional problem must be solved. 

In the present study, a rapid, mixed explicit/implicit 
numerical scheme’ is used to solve the time-dependent thin- 
layer approximation of the three-dimensional Navier-Stokes 
equations for supersonic laminar flow over a hollow cylinder 
with a flare. The use of a hollow cylinder avoids the com- 
plication of prescribing the upstream boundary conditions. 
The objectives of this study are to investigate the three- 
dimensional shock-wave and boundary-layer interaction, to 
examine the associated three-dimensional separated flow 
structure, and to assess the role of circumferential com- 
munication in this interaction flowfield. The experiments 
conducted by Robinson are selected for comparison at 
Mach number 2.8 and unit Reynolds number 0.3168 x 10^ /cm 
with various incident angles. Good agreement is obtained 
between the present computed results and experimental 
measurements of surface pressure and normal force 
distribution. The existence of a pair of vortices is also clearly 
observed. 
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Fig. 3 Fluorescent oil studies on a cylinder with 23 deg cone 
frustrum. 
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Analysis 


Thiii'Layer Approximation 

Figures 4a and 4b show cross-sectional views of the 
computational domain in the (x,r) and {r,<t>) planes. Here, $ 
is the cone half-angle of the flare. The flow is assumed to be 
symmetric in <f> and hence only the half plane, 0 deg 180 
deg, is needed. The mesh is uniformly spaced in both x and <l> 
directions. In the r direction, a fine mesh spacing is used in the 
region near the body, :Sr<ry, to resolve the viscous forces, 
and a coarse mesh spacing is used in the outer region, 
r^<r<r^, where viscous effects are negligible. 

The basic equations of the present analysis are the time- 
dependent compressible Navier-Stokes equations. In high 
Reynolds number flows, however, the viscous effects are 
confined near the wall boundary and are dominated by the 
viscous terms associated with the derivative in the direction 
normal to the wall. The viscous terms with the derivatives 
along the body are comparatively small and insignificant. 
This concept was first discussed by Prandtl in the develop- 
ment of boundary -layer theory and has been applied and 
extended to various techniques, such as the parabolized 
Navier-Stokes equations, or higher order boundary-layer 
theory for many complex flowfield calculations (cf., Refs. 7, 
8, 11-13). The recent development of the thin-layer ap- 
proximation*^ is based on that concept, with retention of all 
the unsteady and nonlinear inviscid terms of the Navier- 
Stokes equations. The retention of time-dependency in the 
governing equations, coupled with the direct inclusion of all 
inviscid terms in the calculation, permits the solution to 
progress naturally from an arbitrary initial condition to an 
asymptotic steady state and, hence, also permits an in- 
teraction flowfield to develop as time progresses. This ap- 
proximation has been tested in the calculation of several, 
flowfield problems and the results confirm its validity and 
applicability.*^'*^ Consequently, the thin-layer approximation 
is used in the present study: All viscous terms associated with 
derivatives along the body are neglected while those with 
second-order derivatives across the thin layer alone are 
retained. Written in weak conservative form in transformed 
cylindrical coordinates, the thin-layer approximated Navier- 
Stokes equations are as follows: 




Fig, 4 Cross sections of the computational domain, a) (x,r) plane; b) 
(r,^) plane. 
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where 


^=x, i7 = r-r^(x), v'^v-mandi 
ian$i=r'i,(x)=0 for x^L 

= tan3 for x>L 


Here {u,v,w) are velocity components in axial, radial, and 
azimuthal, (jr,r,<^), directions; p,/7, and E are density, 
pressure, and total energy, with £’=^, -1-0.5 («^ + ) , 
where is the specific internal energy. The perfect gas 
relations are /?= (7- \)pe^. The flow is assumed laminar and 
the molecular viscosity p is evaluated by Sutherland’s for- 
mula. 
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Time-Split Numerical Technique 

The governing equations, Eqs. (I), are split into three 
groups corresponding to the coordinate directions for 
treatment by three operators, viz., 


drU BrF ^ 


BrU IBrH 

4* 

B( rB<i> 
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where 
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0 


’ 0 ' 

-ptan^j 


pl&ndf 


0 

0 

Qi = 

-p—pw^ 

Qv = 


0 


pvw 



0 


0 


0 


The conventional two-step, explicit MacCormack scheme*^ is 
used in the and operators and in the operator for the 
outer coarse mesh region. In the fine-mesh inner region near 
the body, the operator is split further into hyperbolic, 
and parabolic, operators for the inviscid and viscous 
terms, 


drU 



dt 


(2) 

drU 

drG^ ^ 


dt 

3, 

(3) 


Substitution of (7 and G in Eq. (2) yields 


and is a viscosity parameter, equal to /i for elements 
and equal to (2^ + X) for elements v' and 

and equal to 'ypifPr for element e^. The continuity equation 
has no viscous term and, hence, is trivial to solve. It implies 
that density is treated as a constant in time in the 
operator. To avoid using a costly block-tridiagonal 
procedure, the total energy equation is split into three kinetic 
energy and equations and an internal energy 

equation. The kinetic energy equations are obtained by 
multiplying the corresponding momentum equations by 
u\v \ and w, respectively, and the internal energy equation is 
obtained by subtracting the three kinetic energy equations 
from the total energy equation. Equation (5) is a parabolic 
model equation of diffusion, and can be efficiently solved by 
either a Crank-Nicolson or a first-order implicit difference 
with a simple tridiagonal inversion. The right-hand side of 
Eq. (5) is evaluated explicitly. In the present study, the u\v\ 
and w are solved first and the new values of and w arc 

then used to evaluate the dissipation terms and Aj 

which in general may not be small. Details of the numerical 
method are discussed in Ref. 9. 

The complete numerical procedure is expressed as: 

( A/)I^ ( A/)L, ( 2 At)L^{At)L^ {At) 

Outer region 


/ dv^ 
\ dt 


-\-v' 



-htan^^,) 


di) 



( 



dr) 


) 


+ yp 


dv' 

drj 


ypv' 

r 


(4) 


Note that adding a term ptan^, in Qj and subtracting it in 
is for convenience in the derivation of Eq. (4). 

In the operator, Eq. (4) is solved for v' and p by the 
method of characteristics and then Eq. (2) is integrated in time 
from / to t-k-At by the conventional predictor-corrector 
scheme with v' and p as known values at t + ^ A/. 

To compute the effect of the viscous terms in the inner 
region, Eq. (3) can be rewritten as 


XL,* (A/^) 


Inner region 


>:L^,(Atj)L^^(Atj)L^{ 2 Atj)L^( 2 Aty)L^p{Atj)L^^{Atj) 


Inner region 


xU* 

by which the solution is advanced two time steps from t to 
r-j-2Ar, where 2A/^ = Ar. The maximum stable time step, At, is 
determined by the CFL condition of the explicit operators. 
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Boundary Conditions 

The coordinates are body-oriented and the flow at the 
upstream boundary is assumed to be uniform at supersonic 
freestream conditions ( ) . The angle of attack 
is imposed by setting the upstream velocity components as 
w = w«cosQf,i; = w«sinofCOs<^, and w = w«sinasin0. The 
downstream and far-field boundaries are positioned far 
enough from the three-dimensional interaction region that 
zero-gradient boundary conditions can be used. The sym- 
metry condition is applied at </> = 0 and 180 deg planes. The 
wall is assumed impermeable and no-slip boundary conditions 
are used. The wall is taken to be either isothermal or 
adiabatic, and the wall pressure is evaluated using 

dp 

7 “ =0 at = 0 

dv 

For our present case, the first mesh point is so close to the 
wall that this pressure condition is appropriate. During the 
calculations within the inner mesh, the flux and stress terms at 
the internal boundary, r=rj{x), arc saved, and their net 
quantities are then used as boundary conditions for the outer 
mesh flow calculations. 

Results and Discussion 

The experiments selected for comparison were conducted 
by Robinson.'® The flow conditions arc A/«,=2.8; 
Re I =0.8 X 10^; 15 deg; D= 101 .6 mm; L = 252.5 mm; and 

a = 0, 1, and 4 deg. The wall is assumed adiabatic. A grid of 



AUGUST 1980 


SUPERSONIC LAMINAR FLOW OVER AN INCLINED BODY 


925 


o EXPERIMENTAL DATA (ROBINSON) 
- CALCULATED RESULTS 


INVISCID WEDGE 
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M„ « 2.8, Rbl “ 0.8 X 10®, a - 1® 



M^ ■ 2.8, Re^ - 0.8 X 10®, a • 4® 



Fig. 6 Comparison of normal force distribution, a) a=l deg; b) 
as 4 deg. 


dx//> 


T D JO 


2p 


cos^d0 ) 


(45 X 34 X 20) encloses a computational domain of 


{/L: 

from -0.0156 

to 1.3702 

n/D: 

from 0 

to 1.2338 

<t>: 

from 0 

to 180 deg 



Fig. 5 Comparison of surface pressure, a) a = 0 deg; b) a - 1 deg; c) 
a - 4 deg. 


A geometrically stretched fine mesh, 20 points in ?;-direction, 
is used near the body to t? /D = 0.1397, with a minimum mesh 
spacing (Ajj)„i„ =0.5L/V/?c^. 

Figures 5a*5c show comparisons of the present computed 
results with experimental measurements of surface pressure 
for three incident angles: a = 0, I , and 4 deg. The agreement is 
very good. Also shown in Fig. 5a are the inviscid solutions for 
a sharp cone, two-dimensional wedge, and cylinder-flare, 
respectively. In Fig. 5b, two meridian angles, <^ = 0 and 180 
deg, and in Fig. 5c, three meridian angles 0 = 0, 90, and 180 
deg, are shown. The inviscid solutions shown in Fig. 5c were 
calculated by adding a 30 deg cone in front of the cylinder and 
arc presented to indicate the viscous effects at two meridian 
angles ,0 = 0 and 1 80 deg . 

Note that, at an angle of attack, the “effective** flare angle 
is larger on the windward and smaller on the leeward. If there 





F{g. 8 Unwrapped velocity plots (u ' and w) for a = 4 deg. 


were no circumferential communication or crossflow, the 
streamwise recirculation should be larger on the windward 
and smaller on the leeward. The existence of circumferential 
communication by means of convention and diffusion of 
crossflow, however, makes the Howfield behave oppositely. 
The crossflow sweeps much of the low-momentum fluid from 
the windward boundary layer to the leeward. This process 



results in thinning and strengthening the windward boundary pig. lo Three-dimensional plots of skin friction along streamwise 
layer, while at the same time thickening and weakening the direction. 


leeward boundary layer. Near the cylinder-flare juncture the 


communication is more efficient and effective. While the flare 
shock retards the flow, inducing streamwise recirculation, the 
crossflow is stronger and the swept fluid indeed can be 
carrying “negative” momentum. Overall, it leads to a for- 
ward movement of separation on the leeward and an aft 
movement on the windward, as can be seen in the pressure 
rises and plateaus shown in Figs. 5b and 5c. Similar effects 
were also observed by Deep’’ in an experimental study of 


leeward surface pressure rises gradually and shows a sub- 
stantial difference from the inviscid solution, while the 
windward surface pressure rises sharply and shows less 
viscous effects. This results in a peak of normal force 
distribution near the flare juncture, as shown in Figs. 6a and 
6b. Also shown in the figures arc inviscid solutions using the 
Syvertson and Dennis’® shock expansion method. At small 


shock-induced separation on a cone-cylinder frustum. angle of attack, a = 1 deg, the surface pressure on the leeward 

It is interesting to note that, due to the thickeTiing of the can be higher than that on the windward (Fig. 5b) near the 


boundary layer and the larger separation on the leeward, the streamwise separation point on the lee side, which causes a 
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a) AT (X - L)/q » 0.660 



b) AT (X - D^d = 0.770 

Fig. II Plots of static pressure on (r-^) planes at two x- locations. 




Fig. 12 Plots of static pressure on {x-r) planes at three meridian 
angles. 

small negative normal force distribution in front of the flare 
(Fig. 6a). In general, the present computed results closely 
predict the basic features observed by the experiments. 

Figures 7 and 8 show unwrapped velocity plots, u' vs w, at 
several radial levels for two incident angles, respectively. They 
clearly indicate the variation of streamwise recirculation with 
meridian angle. 

■ One of the most striking features is that, due to the com- 
bination of streamwise recirculation and lee-side crossflow 
separation, there exists a pair of vortices which become 
stretched and are shed into the “wake.” At small incidence, 
oi = I deg, the vortex is weak and not clear. At a = 4 deg, a pair 
of vortices is clearly formed and the feature is very similar to 


the experimental observation of Ericsson et al.’ shown in Fig. 

2. A similar vortex pattern has also been observed on the nose 
of an inclined hemisphere-cylinder by Hsieh and Wang. 

The variations of surface pressure and streamwise skin- 
friction coefficient are shown in Figs. 9 and 10. Note that 
both Figs. 10a and 10b indicate a small secondary recir- 
culation in the tj direction around the lee side. This type of 
secondary recirculation has been observed in many three- 
dimensional experiments; see, for instance. Ref. 20. But in the 
present study, since there is a small pressure oscillation (see 
Figs. 9a and 9b) associated with numerical error due to body 
discontinuity and sharp pressure rise near the juncture, the 
observation is inconclusive and further investigation with 
better mesh resolution is needed. Figure II shows plots of 
static pressure fields on {r-(p) planes at two a: locations and 
Fig. 12 shows those on {x-r) planes at three meridian angles: 

<^ = 0, 90, and 180 deg. At <^ = 90 and 180 deg (Figs. 12b and 
12c) a strong overexpansion, coupled with the leading edge 
and flare-shock compression, is clearly shown. 

The three-dimensional program is coded in such a way as to 
permit treatment of a general body of revolution, with 
r^=f{x). In the present numerical procedure for com- 
putation based on the thin-layer approximation, the program 
requires 0.00048 s per grid point per time-step on a CDC 76(X). 
With a grid of 45x34x20, it takes about 1,25 h of com- 
putation time for a converged steady solution. The option of 
solving the complete Navier-Stokes equations is also coded; it 
requires about 25^o more computation time. For the present 
cases, results from the thin-layer approximation and the 
complete Navier-Stokes equations agree within 3<7o. 

Concluding Remarks 

A supersonic laminar flow over an inclined body of 
revolution has been numerically calculated. The time- 
dependent compressible Navier-Stokes equations, with and 
without a thin-layer approximation, were solved by a rapid, 
mixed explicit-implicit numerical scheme. Good agreement 
was obtained between computed and experimental results for 
surface pressure at different meridian angles and normal 
force distribution. The results show that the crossflow: I) 
effectively transfers low or negative momentum fluid to the 
leeward; 2) causes small separation on the windward and large 
separation on the leeward; 3) leads to a sharp pressure rise on 
the windward, and a smooth turning of streamlines and a 
gradual pressure rise on the leeward; and 4) results in a peak 
in the normal force distribution near the cylinder-flare 
juncture, and a pair of vortices on the leeward. 

References 

* Ericsson, L. E., Reding, P., and Guenther, R. A., “Analytic 
Difficulties in Predicting Dynamic Effects of Separated Flow,” 
Journal of Spacecraft and Rockets^ Vol. 8, Aug. 1971, pp. 872-878. 

^Chyu, W. J. and Coe, C. F., private communication, Moffett 
Field, Calif.. 

^Kuehn, D. M., “Laminar Boundary-Layer Separation Induced by 
Flare on Cylinder at Zero Angle-of- Attack,” NASA TR R-146, 1962. 

^Nielsen, J. H., Lynes, L. L., and Goodwin. F, K., “Theory of 
Laminar Separated Flows on Flared Surfaces Including Supersonic 
Flow with Heating and Cooling,” AGARD Conference Proceedings 
No. 4, Separated Flow, 1966. 

^Lces, L. and Reeves, B. L., “Supersonic Separated and Reat- 
taching Laminar Flow. I, General Theory and Application to 
Adiabatic Boundary-Layer/Shock-Wave Interaction,” AIAA 
Journal, Vol. 2, Nov. 1964, pp. 1907-1920. 

^Stetson, K. F., “Boundary-Layer Separation on Slender Cones at 
Angle-of-Attack,” AIAA Journal, Vol. 10, May 1972, pp. 642-648. 

Lubard, S. C. and HelUwell, W. S., “Calculation of the Flow on a 
Cone at High Angle-of-Attack,” AIAA Journal, Vol. 12, Feb. 1971, 
pp. 277-284. 

®Lin, T C. and Rubin, S. G, “Viscous Flow Over a Cone at 
Moderate Incidence, Pan. 2: Boundary Layer,” Journal of Fluid 
Mechanics, Vol. 59, July 1973, pp. 593-620. 





928 


C. M. HUNG 


AIAA JOURNAL 


’MacCormack, R. W., “An Efficient Numerical Method for 
Solving the Time-Dependent Compressible Navier-Stokes Equations 
at High Reynolds Number,” NASA TM X-73-129, July 1976. 

Robinson, M. L., “Boundary Layer Effects in Supersonic Flow 
over Cylinder-Flare Bodies,” Australian Defense Scientific Service, 
Weapons Research Establishment, Salisbury, South Australia, WRE- 
Repi.-1238 (WR&D), July 1974. 

"Cheng, H. K., Chen. S. Y., Mobley, R., and Huber, C. R., “The 
Viscous Hypersonic Slender-Body Problem: A Numerical Approach 
Based on a System of Composite Equations,” The Rand Corp., Santa 
Monica, Calif., RM 6193-PR, May 1970. 

’^McDonald, H. and Briley, W. R., “Three-Dimensional 
Supersonic Flow of a Viscous or Inviscid Gas,” Journal of Com- 
putational Physics ^ Vol. 19, Oct. 1975, pp. 150-178. 

"Davis, R. T., “Numerical Solution of the Hypersonic Viscous 
Shock-Layer Equations,” AIAA Journal, Vol. 8, May 1970, pp. 843- 
851. 

"Pulliam, T. H. and Steger, J. L., “Implicit Finite-Difference 
Simulation of Three-Dimensional Compressible Flow,” AIAA Paper 
78-10, Huntsville, Ala., Jan. 1978. 


"Baldwin, B. S, and Lomax, H., “Thin Layer Approximation and 
Algebraic Model for Separated Turbulent Flows,” AIAA Paper 78- 
257, Huntsville, Ala., Jan. 1978. 

"MacCormack, R. W. and Baldwin, B. S., “A Numerical Method 
for Solving the Navier-Stokes Equations with Application to Shock- 
Boundary Layer Interactions,” AIAA Paper 75-1, Pasadena, Calif., 
Jan. 1975. 

"Deep, R. A., “Experimental Study of Shock-Induced Separation 
on Conc-Cylindcr-Conical Frustum Body of Revolution,” Redstone 
Arsenal Rept. 6M86, Huntsville, Ala., Dec. 1955. 

"Syvertson, C. A. and Dennis, D. H., “A Second-Order Shock- 
Expansion Method for Applicable to Bodies of Revolution Near Zero 
Lift,” NACA Rept. 1328, 1957. 

*^Hsieh, T. and Wang, K. C., “Concentrated Vortex on the Nose 
of N Inclined Body of Revolution,” AIAA Journal, Vol. 14, May 
1976, pp, 698-700. 

^^Bachalo, W. D., “Experiments on Supersonic Boundary-Layer 
Separation in Three Dimensions,” Journal of Applied Mechanics, 
Transactions of A SME, Series E, Vol. 42, June 1975, pp. 289-294. 



